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Die Erfindung betrifft ein Verfahren zur Zuordnung eines 
Pulsverlaufes zu einem von einer Mehrzahl von Pulstypen 
verschiedener Abklingzeit, insbesondere eines Pulsverlau- 
fes eines Szintillationsdetektors mit zumindest zwei 
Szintillationsmaterialien unterschiedlicher Abklingcha- 
rakteristik. Des weiteren betrifft die Erfindung eine 
Vorrichtung zur Durchf Uhrung des Verfahrens. 

In der Nuklearmedizin werden radioaktiv markierte Sub- 
stanzen zur in-vivo Messung physiologischer Prozesse ver- 
wendet. Eine empf indliche nuklearmedizinische Methode ist 
die Positronen-Emissions-Tomographie (PET). Dabei werden 
isotope verwendet, die unter Aussendung eines Positrons 
2erf alien. Ein besonderer Reiz dieses Verfahrens besteht 
darin, da/3 einige der physiologisch interessantesten Ele- 
mente, wie Kohlenstoff, Stickstoff oder Sauerstoff kurz- 
lebige radioaktive Isotope haben, die beim Zerfall ein 
Positron emittieren. Das ausgesandte Positron annihiliert 
im Gewebe innerhalb einer Strecke von 1-2 mm mit einem 



Elektron und erzeugt zwel, in entgegengesetzte Richtungen 
ausgesandte Gammaquanten einer Energie von 511 keV. Die 
gleichzeitige Detektion der Gammaquanten in gegenuberlie- 
genden Detektoren^ die ublicherweise in einem Detektor- 
ring angeordnet sind^ bildet die Grundlage fur eine tomo- 
graphische Bildrekonstruktion. 

Jeder Detektor des Rings enthait einen Szintillationskri- 
stall/ der beim Auftreffen eines Gammaquants einen Licht- 
blitz mit einer flir das Szintillationsmaterial charakte- 
ristischen Abklingzeit aussendet. Beispielsweise weist 
das oft verwendete Wismutgermanat (BSO) eine Abklingzeit 
von 30 0 ns auf, was zusammen mit der relativ geringen 
Lichtausbeute die erreichbare Ortsauf losung beschrankt. 
Neuere Szintillationsmaterialien mit groBem Potential flir 
PET-Anwendungen umf as sen Lutet ium- Oxy or thos i 1 ikat ( LSO ) 
mit einer 8-fach hoheren Lichtausbeute als BSO und einer 
Abklingzeit von 40 ns, und LuAP (LuAlOg) mit einer Ab- 
klingzeit von 17 ns. Die Lichtpulse der Szintillatorkri- 
stalle werden von einem Photodetektor , typischerweise ei- 
nem Photovervielf acher oder einer Photodiode erfaBt, der 
zur Weiterverarbeitung geeignete elektrische Signale aus- 
gibt. 

Moderne Detektoren enthalten zur Steigerung der Ortsauf- 
15sung zwei oder mehr verschiedene Szintillationskristal- 
le mit unterschiedlichen Abklingzeiten, die mit demselben 
Photodetektor gekoppelt sind. Durch eine Analyse der Aus- 
gabesignale des Photodetektors kann aufgrund der ver- 
schiedenen Abklingzeiten bestimmt werden, welcher Detek- 




torkristall angesprochen hat, und somit die Wechselwir- 
kungstiefe (depth of interaction, DOI) im Detektor erfaJJt 
werden. 



Um diese Zuordnung der Photodetektor-Ausgabesignale zu 
einera der Szintillationskristalle durchzuf uhren, wird iib- 
licherweise ein integrierender Verstarker beim Auslesen 
der Photodetektorsignale eingesetzt. Ein Vergleich der 
Anstiegszeit (definiert als die Zeit von 10% bis 90% des 
Maximalwerts) erlaubt dann eine Untetscheidung der ver- 
schiedenen Detektorpulse. 

in einem anderen Ansatz wird fur jedes Detektorsignal ein 
Startwert (Start des Pulses) und ein Stopwert (Ende des 
Pulses, bzw erster Nulldurchgang) erzeugt. Ein Zeit-zu- 
Amplituden-wandler liefert dann ein Ausgangssignal mit 
einer Amplitude, die der Zeitdauer zwischen dem Start- 
und Stopwert proportional ist, und das zur Diskriminie- 
rung der Detektorpulse verwendet werden kann. 

Allen genannten Unterscheidungs- und Zuordnungsverf ahren 
ist gemein, daB sie eine niedrige ZuverlMssigkeit bei 
verrauschten Signalen aufweisen und eine hohe Abtastrate 
fUr die zuzuordnenden Signalpulse erfordern. 

Der Erfindung, liegt daher die Aufgabe zugrunde, ein Ver- 
f ahren der eingangs genannten Art anzugeben, das eine im 
Hinblick auf Abtastrate effiziente und zuveriassige Zu- 
ordnung eines Pulsverlaufs zu einem von einer Mehrzahl 
von Pulstypen verschiedener Abklingzeit ermoglicht. 



Diese Aufgabe wird erf indungsgemajJ durch ein Verfahren 
zur Zuordnung eines Pulsverlauf es zu einem von einer 
Mehrzahl von Pulstypen verschiedener Abklingzeit, insbe- 
sondere eines Pulsverlauf es eines Szintillationsdetektors 
mit zumindest zwei Szintillationsmaterialien unterschied- 
licher Abklingcharakteristik gelost, das die Verfahren- 
schritte vimfaflt: 

Aufnehmen eines Ausgangspulsverlauf es und Wandeln des 
Pulsverlauf es in ein elektrisches Signal, dessen 
Amplituden-Zeit-Verlauf den Pulsverlauf des Ausgangs- 
pulses reprasentiert ; 

Trans formieren des Amplituden-Zeit-Verlauf s in den 
Frequenzraum um einen den Ausgangspuls reprSsentie- 
renden Amplituden-Frequenz-Verlauf zu erhalten; 

Normieren des Amplituden-Frequenz-Verlauf s um einen 
normierten Amplituden-Frequenz-Verlauf zu erhalten; 

Vergleichen des normierten Amplituden-Frequenz- 
Verlauf s mit einem vorbestimmten Referenzverlauf ; und 

Zuordnen des Ausgangspulsverlauf s zu einem der 
Pulstypen auf Grundlage des Ergebnisses des Ver- 
gleichs . 

Die Erfindung beruht damit auf dem Gedanken, die Zuord- 
nung zu einem der Pulstypen nicht im Zeitraum, sondern im 




Freguenzraum durchzufiihren. Dazu werden die gemessenen 
PulsverlSufe nach ihrer Transformation normiert und mit 
einem Ref erenzverlauf verglichen, der gleichsam einen 
Grenzverlauf zwischen zwei verschiedenen Pulstypen dar- 
stellt, Aus dem Ergebnis dieses Vergleichs kann dann der 
gemessene Pulsverlauf einem der Pulstypen zugeordnet wer- 
den. Dabei erfordert ein Unterscheidung zwischen zwei 
Pulstypen in der Kegel einen Vergleich. Wird eine Zuord- 
nung zu einem von mehr als zwei Pulstypen vorgenommen, 
konnen auch mehrerer Vergleiche mit mehr als einem Refe- 
renzverlauf erf order lich sein. 

Die Ausfuhrung des Vergleichs im Freguenzraum kann bei 
hoherer Zuverlassigkeit wesentlich effektiver (geringere 
Abtastrate) erfolgen als im Zeitraum. Zugleich stehen ef- 
f iziente Algorithmen zur Verfugung um den Zeitverlauf des 
gemessenen Pulses in den Freguenzraum zu trans formier en. 

In einer bevorzugten Ausgestaltung des Verfahrens wird 
das elektrische Signal vor der Transformation in den Fre- 
guenzraum einer Analog/Digital-Wandlung mit einer vorbe- 
stimmten Abtastrate unterzogen, um einen den Ausgangspuls 
reprasentierenden^ diskreten Amplituden-Zeit-Verlauf zu 
erhalten. 

Zweckmaflig wird bei dem erf indungsgemSBen Verfahren bei 
dem Schritt des Trans formierens der diskrete Amplituden- 
Zeit-Verlauf einer diskreten Fouriertrans formation (DFT) 
unterzogen, um einen diskreten Amplituden-Freguenz- 
Verlauf zu erhalten. Insbesondere kann die diskrete Fou- 
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riertrans formation mit einer schnellen Four iertrans forma- 
tion (FFT) berechnet werden. 

Bei dem Schritt des Normierens wird vorteilhafterweise 
der Amplituden-Freguenz-Verlauf auf die Amplitude bei der 
Frequenz Null bezogen, insbesondere wird der Amplituden- 
Frequenz-verlauf zweckmSBig durch die Amplitude bei der 
Frequenz Null geteilt. 

in einer bevorzugten Weiterbildung des erf indungsgemSBen 
Verfahrens wird bei dem Schritt des Vergleichens der Dif- 
ferenz-verlauf des normierten Amplituden-Frequenz- 
Verlaufs und des Ref erenzverlauf s gebildet, der Diffe- 
renz-verlauf fUr jede Frequenz mit einem vorbestimmten 
Wichtungsfaktor multipliziert, und als Zuordnungsparame- 
ter die Summe dieser Produkte iiber alle Frequenzen be- 
rechnet . 



Dabei kann die Zuordnung zu einem der Pulstypen auf 
Grundlage des Vorzeichens und/oder des Betrags des Zuord- 
nungsparameters erfolgen. Insbesondere kann im Fall, wo 
der Ausgangspulsverlauf einem von zwei Pulstypen ver- 
schiedener Abklingzeit zugeordnet wird, die Zuordnung auf 
Grundlage nur des Vorzeichens des Zuordnungsparameters 
erfolgen. Dies gestattet einen besonderes einfache Ent- 
scheidung Uber die Zuordnung eines gemessene Pulses zu 
einer der Pulstypen. 

Der Wichtungsfaktor ist in einer Ausgestaltung des Ver- 
fahrens der art gewShlt, dal3 er ein Maximum bei kleinen 



Frequenzen hat und zu groJJen Frequenzen hin abnimmt. Dies 
tragt der Beobachtung Rechnung, dafl die niedrigeren Fre- 
quenzen starker als hohere Frequenzen zur Diskriminierung 
der Pulstypen beitragen. 

Dabei ist es im vorstehend erlauterten Zusammenhang be- 
vorzugt, wenn die hochf requenten Komponenten des elektri- 
schen Signals vor der Transformation in den Frequenzraum 
durch ein TiefpaBfilter unterdrUckt werden. Hierdurch 
kdnnen stSrende RUckf altungen von Signalen, deren Fre- 
quenz hSher als die halbe Abtastrate liegt, bei der Digi- 
talisierung (ADC) des Signals vermieden werden. 

Das Verfahren kann in besonders vorteilhafter Weise da- 
durch weitergebildet sein, da/3 vor dem Aufnehmen zuzuord- 
nender PulsverlSufe ein Ref erenzverlauf bestimmt wird 
durch 

- Aufnehmen von EichpulsverlSufen mit bekannter Zuord- 
nung zu zwei Pulstypen mit verschiedener Abklingzeit 
und Wandeln jeweils des Pulsverlauf es in ein elektri- 
sches Signal, dessen Amplituden-zeit-Verlauf den 
Pulsverlauf des Eichpulses reprSsentiert; 

- Transformieren jeweils des Amplituden-Zeit-Verlauf s 
in den Frequenzraum um einen den Eichpuls repr^sen- 
tierenden Araplituden-Frequenz-Verlauf zu erhalten; 



Normieren des Amplituden-Frequenz-Verlauf s um einen 
normierten Amplituden-Frequenz-Verlauf zu erhalten; 
und 

Festlegen eines Ref erenzwertes fur jede Frequenz der- 
art, daB die Amplitudenwerte der Eichpulse des ersten 
Puistyps fur diese Frequenz im wesentlichen groBer 
als der Ref erenzwert, und die Amplitudenwerte der 
Eichpulse des zweiten Puistyps flir diese Frequenz im 
wesentlichen kleiner als der Referenzwert sind. 

Insbesondere kann der Referenzwert fUr jede Frequenz mit 
Vorteil dadurch festgelegt werden, daB jeweils ein Ampli- 
tuden-Histogramm fur die verschiedenen Pulstypen aufge- 
stellt wird, der Schnittpunkt der Einhullenden der Histo- 
gramme bestimmt wird, und der bestimmte Schnittpunkt als 
der Referenzwert festgelegt wird. 

Als besonders vorteilhaft wird es im vorstehend erlSuter- 
ten Zusammenhang angesehen, wenn vor dem Aufnehmen zuzu- 
ordnender Pulsverlaufe der Wichtungsf aktor fur jede Fre- 
quenz bestimmt wird durch 

Bestimmen des mittleren Betrags der Abweichung der 
Amplitudenwerte von dem Referenzwert fur diese Fre- 
quenz flir jeden der Pulstypen; und 

Festlegen eine hohen Wichtung flir diese Frequenz bei 
einem groBen mittleren Abweichungsbetrag und einer 
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niedrigen Wichtung fUr diese Frequenz bei einem nied- 
rigen mittleren Abweichungsbetrag. 

Eine erf indungsgeinaJ3e Vorrichtung zur Durchfuhrung des 
beschriebenen Verfahrens umfaBt 

ein Mittel zum Aufnehitien eines Ausgangspulsverlauf es 
und einem Mittel zum Wandeln des Pulsverlauf es in ein 
elektrisches Signal , dessen Amplituden-Zeit-Ver lauf 
den Pulsverlauf des Ausgangspulses reprasentiert; 

ein Mittel zum Trans formier en des Amplituden-Zeit- 
Verlaufs in den Prequenzraum urn einen den Ausgangs- 
puis reprasentierenden Amplituden-Frequenz-Verlauf zu 
erhalten ; 

ein Mittel zum Normieren des Amplituden-Frequenz- 
Verlaufs urn einen normierten Amplituden-Frequenz- 
Verlauf zu erhalten; 

ein Mittel zum Vergleichen des normierten Amplituden- 
Frequenz-verlaufs mit einem vorbestimmten Referenz- 
verlauf ; 

ein Mittel zum Zuordnen des Ausgangspulsverlauf s zu 
einem der Pulstypen auf Grundlage des Ergebnisses des 
Vergleichs und zur Ausgabe des Zuordnungsergebnisses; 

ein Mittel zum Bestimmen des Ref erenzverlauf es; 
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ein Mittel zum Bestimmen der Wichtungsf aktoren. 

Dabei konnen die Mittel grundsStzlich jeweils sowohl 
durch eine Hardwarekomponente , wie einen anwenderprogram- 
mierbaren Logikschaltkreis (FPGA), durch Software in ei- 
nem geeignet eingerichteten Rechner, oder durch Firmware 
gebildet sein. Ebenso konnen die Mittel in einem digita- 
len Signalprozessor (DSP) implementiert sein. 

Weitere vorteilhafte Ausgestaltungen, Merkmale und De- 
tails der Erfindung ergeben sich aus den abhSngigen An- 
spruchen, der Beschreibung des Ausfuhrungsbeispiels und 
den Zeichnungen. 

Die Erfindung soli nachfolgend anhand eines Ausfuhrungs- 
beispiels im Zusammenhang mit den Zeichnungen naher er- 
lautert werden. Es sind jeweils nur die fur das VerstSnd- 
nis der Erfindung wesentlichen Elemente dargestellt. Da- 
bei zeigt 

Figur 1 eine Blockdiagramm einer Vorrichtung zur Zuord- 
nung eines gemessenen Pulsverlauf es nach einem 
Ausflihrungsbeispiel der Erfindung; 

Figur 2 eine Auftragung von Amplituden-Zeit-Verlauf en 
flir je zwei Signalpulse von zwei verschiedenen 
Szint illationskris tallen ; 

Figur 3 eine Auftragung der Amplituden-Frequenz- 
Verlaufe der Signalpulse von Fig. 2; 



eine Auftragung der normierten Amplituden- 
Frequenz-VerlSuf e der Signalpulse von Fig. 2 
und eines Ref erenzverlauf s ; 

eine Auftragung der Dif f erenz . zwischen den nor- 
mierten Amplituden-Frequenz-Verlauf en und dem 
Ref erenzverlauf von Fig. 4; 

ein Diagramm, das die Wichtungsf aktoren fiir die 
Frequenzkomponenten von Fig. 5 angibt; 

eine Auftragung der gewichteten Dif f erenzver- 
lauf e von Fig. 5; und 

eine Auftragung des nach dem erf indungsgem^Ben 
Verfahren bestimmten Zuordnungsparameters fiir 
100 Ausgangspulse von zwei verschiedenen Szin- 
t illationskristallen • 

Ein Beispiel einer Anwendung des erf indungsgemaBen Ver- 
fahren zur Zuordnung eines Pulsverlauf es zu zwei Pulsty- 
pen wird nachfolgend mit Bezug auf die Figuren 1 bis 8 im 
Detail beschreiben. 

Der verwendete Szintillationsdetektor 10 enthalt zwei 
verschiedene Szintillationskristall-Schichten 12 und 14, 
die hintereinander vor einem Photodetektor 16 angeordnet 
Bind. Der Szintillationskristall 12 besteht im Ausfuh- 
rungsbeispiel aus LSO (Lutetium-Oxyorthosilikat ) mit ei- 



Figur 4 



Figur 5 



Figur 6 



Figur 7 



Figur 8 
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ner Abklingzeit von 40 ns^ und der Szintillationskristall 
14 aus LuYAP (LUi.^Y^lOa) mit einer Abklingzeit von 17 ns. 

Zur Bestimmung des Ref erenzverlauf es und der Wichtungs- 
faktoren wurden zunachst Eichpulse mit bekannter Zuord- 
nung zu dem schnellen LuYAP-Kristall oder dem langsameren 
LSO-Kristall auf genommen. Die vom Photodetektor 16 gelie- 
ferten elektrischen Pulsveriaufe der Eichpulse wurden mit 
einem Butterworth-Filter 22 zweiter Ordnung und einer 
Grenzfrequenz von 20 MHz tiefpaB-gef iltert . Dann wurden 
die Pulse mit einer Abtastrate fg von 40 MHz, entspre- 
chend dem doppelten der Grenzfrequenz des TiefpaBf ilters 
von dem A-D-Wandler 24 digitalisiert . 

Das Ergebnis nach Tief paB-Filterung und Digitalisierung 
ist in Fig. 2 fur jeweils zwei Pulse des LSO-Kristalls 12 
(Bezugszeichen 42H, 42N) und des LuYAP-Kristalls 14 
(Bezugszeichen 44H, 44N) gezeigt. Dabei ist zur besseren 
Illustration des Verfahrens jeweils ein Puis mit hoher 
Energie (LSO: 42H, LUYAP: 44H) und ein Puis mit niedriger 
Energie (LSO: 42N, LUYAP: 44N) ausgewahlt und darge- 
stellt . 

Dann werden die Amplituden-Zeit-Veriauf e d.er Zuordnungs- 
vorrichtung 26 eingegeben, die in anwenderprogrammierba- 
ren Logikschaltkreise 28, 30, 32, 34 und 36 umfafit, deren 
Funktion aus der nachf olgenden Beschreibung deutlich 
wird. 



Die Amplituden-Zeit-Verlaufe der Eichpulse wurden zu- 
nachst mittels einer 32-Punkt-FFT f ouriertransformiert 
(Bezugszeichen 28), wodurch bei der genannten Abtastrate 
eine Auflosung Af von 4 0 MHz/32 =1,25 MHz erzielt wird. 

Da sich die Pulsformen der beiden Detektoren vor allem im 
Bereich niedrigerer Frequenzen unterhalb von 10 MHz un- 
terscheiden, wurden auch bei einer niedrigeren Abtastrate 
von fg = 20 MHz, einem TiefpaBf ilter mit einer Grenzfre- 
quenz von 10 MHz und einer 16-Punki:-FFT, bei derselben 
Auflosung Af = 1,25 MHz, gute Zuordnungsergebnisse er- 
zielt werden. Wie weiter unten im Detail beschrieben, er- 
gab der Rest mit der hoheren Abtastrate eine sicherere 
Zuordnung der Pulse von 97% erfordert aber wegen der ho- 
heren Punktzahl der FFT hoheren Aufwand fur die Berech- 
nung. 

Das Ergebnis der 32-Punkt FFT fur die vier Pulsverl'aufe 
von Fig. 2 ist in Fig. 3 dargestellt. Die i-te FFT- 
Komponente auf der horizontalen Achse entspricht dabei fi 
= (i-l)*l,25 MHz, die vertikale Achse gibt die Amplitude 
A(fi) der Pulsveriaufe bei der Frequenz fi an. Wegen der 
Symmetrie der 32 FFT-Komponenten sind in der Fig. 3 nur 
die ersten 16 FFT-Komponenten gezeigt. Die Amplituden- 
Frequenz-VerlSufe der vier Eichpulse sind durch Quadrate 
52H (LSO, hohe Energie), Kreise 54H (LuYAP, hohe Ener- 
gie), Sterne 52N (LSO, niedrige Energie) bzw. Kreuze 54N 
(LuYAP, niedrige Energie) dargestellt. Die unterschiedli- 
che Gesamtenergie ist an den unterschiedlichen DC- 
Amplituden bei der Komponente i=l deutlich zu erkennen. 
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Zur Normierung wird der Amplituden-Frequenzverlauf jeden 
Pulses durch seine erste Frequenzkomponente A{f^) geteilt 
(Bezugszeichen 30). Die resultierenden normierten Ampli- 
tuden-Frequenz-Verlaufe A„(fi) sind dann nicht mehr von 
der Energie der Pulse abhangig. 

Fig. 4 zeig-t die normierten Verlaufe fUr LSO mit Bezugs- 
zeichen 62H und 62N, fur LuYAP mit Bezugszeichen 64H und 
64N. Wie aus der Fig, 4 deutlich zu erkennen, sind die 
Signale der verschiedenen Detektoren im Frequenzraum 
durch eine Liicke voneinander getrennt, die bei den nied- 
rigeren Frequenzkomponenten bis fg besonders ausge- 

pragt ist. Bin Ref erenzverlauf 66, der in der Mitte der 
Liicke zwischen den beiden Pulstypen verlSuft, ist eben- 
falls in Fig. 4 dargestellt. Wie die Figur ebenfalls 
zeigt^ kann mehr als eine Frequenzkomponente flir die Un- 
terscheidung der beiden Pulstypen veirwenden werden. 

Der Ref erenzverlauf 66 wurde fur jede Frequenz f^ dadurch 
erhalten, dai3 die Amplitudenwerte fUr eine Reihe von 
Eichpulsen jeden Pulstyps in ein Amplituden-Histogramm 
eingetragen wurden. Der Schnittpunkt der Einhiillenden der 
jeweiligen Amplituden-Histogramme ergibt den Wert des Re- 
f erenzverlauf s R(fi) fur diese Frequenz f^. 

Nun wurde die Differenz zwischen den normierten Amplitu- 
den-Frequenz-Veriauf en A^, ( f ^ ) und dem Ref erenzverlauf 
R(fi) berechnet (Bezugszeichen 32). Das Ergebnis zeigt 
Fig. 5 fur die LSO-Pulse als 72H und 72N, und fur die Lu- 
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YAP-Pulse als 74H und 74N. Die Dif f erenzwerte sind dabei 
fur die LSO-Pulse grofltenteils positiv, fur die Lu YAP- 
Pulse im wesent lichen negativ. wie aus Fig. 5 ersicht- 
lich, lassen sich die verschiedenen Pulse bei niedrigeren 
Frequenzkomponenten gut trennen, flir hohere Frequenzkom- 
ponenten geht die Unterscheidungskraf t jedoch gegen Null, 

Um diesem Umstand Rechnung zu tragen, wurde flir jede Fre- 
quenzkomponente ein Wichtungsf aktor W(fi) festgelegt, der 
ein MaJ3 fur die Eignung dieser Frequenzkomponente zur 
Pulstypen-Diskriminierung darstellt. Dabei wird jeweils 
der mittlere Betrag der Abweichung der Amplitudenwerte 
von dem Referenzwert flir jeden der Pulstypfen als MaJ3 flir 
die Unterscheidungseignung genommen. Aus den in Fig. 5 
dargestellten Amplituden lassen sich etwa die in Fig. 6 
gezeigten Wichtungsf aktoren W(fi) flir jede Frequenzkompo- 
nente fi bilden. Wie auch anschaulich aus Fig. 5 klar, 
erhalten dabei die Komponenten mit niedrigem Index, ins- 
besondere die Frequenzkomponenten bei 2,5 MHz und 3,75 
MHz, ein hohes Gewicht, wahrend die Wichtungsf aktoren flir 
die nicht zur Diskriminierung beitragenden hoheren Fre- 
quenzen auf Null gesetzt sind. 

Der so ermittelte Ref erenzverlauf R(fi) und die Wich- 
tungsf aktoren W(fi) werden gespeichert, und konnen dann 
flir die Zuordnung gemessener Pulsverlaufe verwendet wer- 
den, deren Ur sprung von einem der beiden Detektoren nicht 
bekannt ist, wie nachfolgend dargestellt. 



Zur Illustration wird die Zuordnung anhand der bisher 
verwendeten Eichpulsen des LSO und des LuYAP Kristalls 
erlautert. Wird ein neuer Puis unbekannter Zuordnung ge- 
messenen, so werden die oben beschriebenen Schritte mit 
Ausnahme der Berechnung des Refer enzverlaufs und der 
Wichtungsf aktoren durchgef lihrt . 

bann wird die Differenz zwischen dem normierten Amplitu- 
den-Frequenz-Verlauf A„(fi) und dem Ref erenzverlauf R(fi) 
fur jede Frequenz mit dem Wichtungsf akt or multipliziert 

D„(fi) = W(fi) * (A«(fi) - R(fi)) 

urn eine gewichtete Differenz Dw(fi) zu erhalten 
(Bezugszeichen 34). Die so gewichteten Werteverlauf e sind 
in der Fig. 7 fUr die LSO-Pulse als B2H bzw. 82N und flir 
die LuYAP-Pulse als 84H und 84N dargestellt. 

Ein Aufsummieren liber die Beitrage von alien Frequenzen 
ergibt dann den Zuordnungsparameter Z fur jeden Puis 
(Bezugszeichen 34). Zu Zuordnung geniigt das Vorzeichen 
von Z. Ein positiver Wert von Z f lihrt zu einer Einordnung 
des gemessenen Pulses als Signal des LSO-Kristalls, ein 
negativer Wert zu einer Einordnung als Signal des LuYAP- 
Kristalls (Bezugszeichen 36). 

Zur Untersuchung der ZuverlSssigkeit des erf indungsgemS- 
Ben Verfahrens wurden nach einmaliger Festlegung des Re- 
f erenzverlauf s R(fi) und der Wichtungsf aktoren W(fi) ein- 
hundert Ausgangspulse, flinfzig des LSO-Kristalls und 
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fiinfzig des LuYAP-Kristalls nach dem vorstehenden erlau- 
terten Verfahren zugeordnet. Figur 8 zeigt jeweils den 
Zuordnungsparameter z fur die LSO-Pulse (Bezugszeichen 
92) und die LuYAP-Pulse (Bezugszeichen 94). Die Z-Werte 
fur die vier Eichpulse sind als 92H, 92N, 94H und 94N 
ebenfalls eingezeichnet. Wie aus der Figur ersichtlich, 
fuhrt das erf indungsgemSBe Verfahren bei 97% der Pulse zu 
einer korrekten Zuordnung. 



AnsprlichQs 

Forsch unaszentrum Julich GmbH. Wilhelm-Johnen-Str . , 
D-52425 Julich 

Verfahren zur Zuordnuna eines Pulsverlauf es zu einem von 
einer Mehrza hl von Pulstypen verschiedener Abklinazeit: 
und Vorrichiiun a zu dessen Durchfuhruna r Pulsf ormdiskri- 
minieruna) 

1. Verfahren zur Zuordnung eines Pulsverlauf es zu einem 
von einer Mehrzahl von Pulstypen verschiedener Ab- 
klingzeit, insbesondere eines Pulsverlauf es eines 
Szintillationsdetektors mit zumindest zwei Szintil- 
lationsmaterialien unterschiedlicher Abklingcharak- 
teristik^ umfassend die Verf ahrenschritte: 

- Aufnehmen eines Ausgangspulsverlauf es und Wandeln 
des Pulsverlauf es in ein elektrisches Signal, des- 
sen Amplituden-Zeit-Verlauf den Pulsverlauf des 
Ausgangspulses reprSsentiert; 

- Trans formier en des Amplituden-Zeit-Verlauf s in den 
Frequenzraum urn einen den Ausgangspuls repr^sen- 
tierenden Amplituden-Frequenz-Verlauf zu erhalten; 

- Normieren des Amplituden-Frequenz-Verlauf s um ei- 
nen normierten Amplituden-Frequenz-Verlauf zu er- 
halten; 




- Vergleichen des normierten Amplituden-Frequenz- 
Verlaufs xnit einem vorbestimmten Ref erenzverlauf ; 
und 

Zuordnen des Ausgangspulsverlauf s zu einem der 
Pulstypen auf Grundlage des Ergebnisses des Ver- 
gleichs . 

2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, 
daB das elektrische Signal vor der Transformation in 
den Frequenzraum einer Analog/Digital-Wandlung mit 
einer vorbestimmten Abtastrate unterzogen wird^ urn 
einen den Ausgangspuls reprasentierenden, diskreten 
Amplituden-Zeit-Verlauf zu erhalten. 



3. Verfahren nach Anspruch 2, dadurch gekennzeichnet , 
daJ3 bei dem Schritt des Trans formier ens der diskrete 
Amplituden-Zeit-Verlauf einer diskreten Fourier- 
transformation (DPT) unterzogen wird^ um einen dis- 
kreten Amplituden-Frequenz-Verlauf zu erhalten. 

4. Verfahren nach Anspruch 3, dadurch gekennzeichnet, 
daJ3 die diskrete Fouriertransf ormation mit einer 
schnellen Fouriertransf ormation (FPT) berechnet 
wird. 



5. Verfahren nach einem der vorigen Anspriiche, dadurch 
gekennzeichnet, dafl bei dem Schritt des Normierens 
der Amplituden-Frequenz-Verlauf auf die Amplitude 
bei der Frequenz Null bezogen wird, insbesondere daB 



der Amplituden-Frequenz-Verlauf durch die Amplitude 
bei der Frequenz Null geteilt wird. 

6. Verfahren nach einem der vorigen Anspriiche/ dadurch 
gekennzeichnet, daJ3 bei dem Schritt des Vergleichens 
der Dif f erenz-Verlauf des normierten Amplituden- 
Frequenz-Verlauf s und des Ref erenzverlauf s gebildet 
wird^ der Dif f erenz-Verlauf fiir jede Frequenz mit 
einem vorbestiimnten Wichtungsf aktor multipliziert 
wird, und als Zuordnungsparameter die Summe dieser 
Produkte uber alle Frequenzen berechnet wird. 

7. Verfahren nach Anspruch 6, dadurch gekennzeichnet , 
daJ3 die Zuordnung zu einem der Pulstypen auf Grund- 
lage des Vorzeichens und/oder des Betrags des Zuord- 
nungsparameters erf olgt . 

8. Verfahren nach Anspruch 6 oder 1, dadurch gekenn- 
zeichnet, dai3 der Ausgangspulsverlauf zu einem von 
zwei Pulstypen verschiedener Abklingzeit zugeordnet 
wird, und die Zuordnung auf Grundlage nur des Vor- 
zeichens des Zuordnungsparameters erf olgt • 

9. Verfahren nach einem der Ansprliche 6 bis 9, dadurch 
gekennzelchnet, daJ3 der Wichtungsf aktor ein Maximum 
bei kleinen Frequenzen hat und zu groflen Frequenzen 
hin abniirant* 

10. Verfahren nach einem der vorigen Ansprliche, dadurch 
gekennzeichnet, daJ3 hochf requente Komponenten des 
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elektrischen Signals vor der Transformation in den 
Frequenzraum durch ein TiefpaBf ilter unterdrUckt 
warden • 

11 • Verfahren nach einem der vorigen Ansprliche, dadurch 
gekennzeichnet, daB vor dem Aufnehmen zuzuordnender 
Pulsveriaufe ein Ref erenzverlauf bestimmt wird durch 

- Aufnehmen von Eichpulsverlauf en mit bekannter Zu- 
ordnung zu zwei Pulstypen mit verschiedener Ab- 
klingzeit und Wandeln jeweils des Pulsverlauf es in 
ein elektrisches Signal, dessen Amplituden-Zeit- 
Verlauf den Pulsverlauf des Eichpulses reprasen- 
tiert; 

- Transformieren jeweils des Amplituden-Zeit- 
Verlaufs in den Frequenzraum um einen den Eichpuls 
reprSsentierenden Amplituden-Frequenz-Verlauf zu 
erhalten; 

- Normieren des Amplituden-Frequenz-Verlauf s um ei- 
nen normierten Amplituden-Frequenz-Verlauf zu er- 
halten; und 

- Festlegen eines Ref erenzwertes flir jede Frequenz 
der art, dafl die Amplitudenwerte der Eichpulse des 
ersten Pulstyps flir diese Frequenz im wesentlichen 
groBer als der Ref erenzwert , und die Amplituden- 
werte der Eichpulse des zweiten Pulstyps fur diese 



Frequenz im wesentlichen kleiner als der Referenz- 
wert sind. 

12. Verfahren nach Anspruch 11^ dadurch gekennzeichnet, 
daJ3 der Referenzwert fur jede Frequenz dadurch fest- 
gelegt wird, daJ3 jeweils ein Amplituden-Histogramm 
fUr die verschiedenen Pulstypen aufgestellt wird, 
der Schnittpunkt der Einhlillenden der Histogranune 
bestimmt wird, und der bestimmte Schnittpunkt als 
der Referenzwert festgelegt wird. 

13. Verfahren nach Anspruch 11 oder 12, dadurch gekenn- 
zeichnet^ daB vor dem Aufnehmen zuzuordnender Puls- 
verlaufe der Wichtungsf aktor fUr jede Frequenz be- 
stimmt wird durch 

- Bestimmen des mittleren Betrags der Abweichung der 
Amplitudenwerte von dem Referenzwert fur diese 
Frequenz flir jeden der Pulstypen; und 

- Festlegen eine hohen Wichtung fiir diese Frequenz 
bei einem groJJen mittleren Abweichungsbetrag und 
einer niedrigen Wichtung fiir diese Frequenz bei 
einem niedrigen mittleren Abweichungsbetrag. 

14. Vorrichtung zur Durchfiihrung des Verfahrens nach ei- 
nem der Anspriiche 1 bis 13, mit 

- einem Mittel zum Aufnehmen eines Ausgangspulsver- 
laufes und einem Mittel (10) zum Wandeln des Puis- 
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verlaufes in ein elektrisches Signal, dessen 
Amplituden-Zeit-Verlauf den Pulsverlauf des Aus- 
gangspulses reprasentiert ; 

- einem Mittel (28) zum Trans formier en des Amplitu- 
den-Zeit-Verlauf s in den Frequenzraum um einen den 
Ausgangspuls reprasentierenden Amplituden- 
Frequenz-Verlauf zu erhalten; 

- einem Mittel (30) zum Normieren des Amplituden- 
Frequenz-Verlauf s um einen normierten Amplituden- 
Frequenz-Verlauf zu erhalten; 

- einem Mittel (32,34) zum Vergleichen des normier- 
ten Amplituden-Frequenz-Verlauf s mit einem vorbe- 
stimmten Ref erenzver lauf ; und 

- einem Mittel (36) zum Zuordnen des Ausgangspuls- 
verlaufs zu einem der Pulstypen auf Grundlage des 
Ergebnisses des Vergleichs und zur Ausgabe des Zu- 
ordnungsergebnisses . 

15. Vorrichtung nach Anspruch 14, bei dem das Mittel 
(32,34) zvim Vergleichen und das Mittel (36) zum Zu- 
ordnen durch einen anwenderprogrammierbaren Logik- 
schaltkreis (FPGA) oder einen digitalen Signalpro- 
zessor (DSP) oder (PC) gebildet ist. 
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16. Vorrichtung nach Anspruch 15, bei dem das Mittel 
(28) zum Trans formieren und das Mittel (30) zum Nor- 
mieren durch einen anwenderprograitunierbaren Logik- 
schaltkreis (FPGA) oder einen digitalen Signalpro- 
zessor (DSP) oder (PC) gebildet ist. 
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Zusanunenf assuna; 

Forschunaszentrum Jlilich GmbH. Wilhelm- Johnen-S-br > , 
D-52425 Julich 

Verfahren zur Zuordnuna eines Pulsverlauf es zu einem von 
einer Mehrzahl von Puls-hypen verschiedener Abklinazei-t 
und Vorrichiiuna zu des sen DurchfUhruna (Pulsf ormdiskri- 
minieruna^ 

Ein Verfahren zur Zuordnung eines Pulsverlauf es zu einem 
von einer Mehrzahl von Pulstypen verschiedener Abkling- 
zeit, insbesondere eines Pulsverlauf es eines Szintillati- 
onsdetektors mit zumindest zwei Szintillationsmaterialien 
unterschiedlicher Abklingcharakteristik, umfaBt die Ver- 
f ahrenschritte : 

- Aufnehmen eines Ausgangspulsverlauf es und Wandeln 
des Pulsverlauf es in ein elektrisches Signal, dessen 
Amplituden-Zeit-Verlauf den Pulsverlauf des Aus- 
gangspulses reprSsentiert; 

Trans formier en des Amplituden-Zeit-Verlauf s in den 
Frequenzraum um einen den Ausgangspuls reprasentie- 
renden Amplituden-Frequenz-Verlauf zu erhalten; 

Normieren des Amplituden-Frequenz-Verlauf s um einen 
normierten Amplituden-Frequenz-Verlauf zu erhalten; 




- Vergleichen des normierten Amplituden-Frequenz- 
Verlaufs mit einem vorbestimmten Ref erenzverlauf ; 
und 

Zuordnen des Ausgangspulsverlauf s zu einem der 
Pulstypen auf Grundlage des Ergebnlsses des Ver- 
gleichs . 

(Figur 1) 
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